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PROPAGATING  MOOES  ON  A BURIED  LEAKY  COAXIAL  CARLE 


by 

Steven  W.  Plate,  David  C.  Chang  and  Edward  F.  Kuester 

1 . Introduction 

The  problem  of  wave  propagation  along  a buried,  insultated  wire  has 
long  been  associated  with  long  distance  communication  at  very  low  and 
extremely  low  frequencies  (VLF  and  OLF).^  Operating  at  Ion/'  wavelengths, 
the  wire  is  usually  very  close  to  the  earth  surface,  say  within  one  or  two 
skin-depths  so  that  the  interaction  between  the  wire  and  the  air- inter  face 
plays  an  important  role  in  determining  the  propagation  constant  of  the  current 

r 2 1 

waves  supported  by  the  wire  structure.  1 More  recently,  buried  wires 
operating  at  the  higher  end  of  the  radio  spectrum,  100  MHz  for  instance, 
have  also  found  important  applications  as  wave  guiding  structures  in  the 
design  of  groundwave  radar  detection  systems  either  for  vehicle  monitoring 
or  perimeter  surveillance.  These  wires  again  have  to  be  placed  close  to  the 
earth  surface  in  order  to  avoid  unnecessary  loss  into  the  surrounding  earth 
ns  well  as  to  enhance  area  coverage.  In  many  practical  applications,  leaky 
coaxial  cables,  either  in  the  form  of  braided  wires  or  periodic  slots  in 
the  outer  sheath,  are  used  to  allow  electromagnetic  waves  to  leak  out 
continuously  from  the  coaxial  region  into  the  surrounding  medium.  It  is 
apparent  that  the  characteristics  of  the  propagating  modes  supported  by 
such  a waveguiding  system  would  Have  to  be  influenced  by  the  air-earth 
interface.  The  purpose  of  this  report  is  to  investigate  such  an  influence. 
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A related  problem  concerning  an  elevated  wire  structure  above  earth 
surface  has  been  studied  extensively.  It  is  found  that  a bare 

or  a dielectrically-coated  wire  usually  supports  two  distinct  modes. 

One  of  the  two  modes  has  a field  structure  more  concentrated  between  the 
wire  and  its  image  and  approaches  the  conventional  transmission- line  mode  in 
the  limiting  case  of  a perfectly-conducting  earth.  Such  a mode  is  referred 
to  as  structure-attached.  The  second  mode,  however,  is  more  spead  out  over 
the  air-earth  interface,  much  like  a groundwnve  field  guided  along  the 
direction  of  the  wire.  This  mode  is  then  referred  to  as  surface-attached. 
Although  the  two  generally  exhibit  very  different  properties,  both  propa- 
gation constants  are  very  close  the  wavenumber  in  air,  and  degeneracy  of 

r 121 

the  two  can  occur  when  the  wire  parameters  are  properly  chosen1  . The 
result  of  our  investigation  indicates  that  such  a phenomenon  does  not 
exist  in  the  case  of  a buried  wire,  however. 


2.  The  Modal  Kquat ion 

Consider  an  infinitely  long  thin  cable  of  exterior  radius  a buried 

in  the  earth  at  a depth  h parallel  to  the  z-axis.  The  earth, region  1, 

(x  > 01  is  assumed  to  be  nonmagnetic,  having  the  refractive  index 

n^  * (rrl  ► ici^/ioc^l  and  region  2,  the  refractive  index  n,.  In  this 

problem,  region  2 usually  represents  free  space  which  means  that  n7 

reduces  to  unity;  however,  n,  is  left  arbitrary  for  generality.  All 

«• 

field  quantities  nre  assumed  to  vary  as  exp (iak^z  -iwt)  where  a is  the 
yet  undetermined,  complex  propagation  constant  of  a discrete  mode.  The 
geometry  of  this  problem  is  illustrated  in  figure  1.  We  further  assume 
that  the  cable  is  thin  compared  to  the  depth  at  which  it  is  buried  (a«h) 
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L. 


and  to  the  wavelength  in  region  1 (k  a|n,|  «1  where  k * u>(u  c ) ). 

o 1 o o o 

This  implies  that  the  current  on  the  cable  will  be  virtually  only  x-directed 
and  will  be  ^-independent . The  corresponding  modal  equation  for  propagating 
modes  supported  by  the  wire  structure  is  then  found  by  assuming  a current  I 
on  the  cable,  calculating  the  average  z -directed  E-field  on  the  exterior 
surface  of  the  cable  due  to  the  current  1,  and  then  enforcing  the  boundary 
condition  Ez  * 12(a)  where  2(a)  is  the  surface  impedance  defined  at  the 
surface  of  the  cable.  All  of  the  dependency  on  the  internal  structure  of 
the  cable  is  contained  in  the  surface  impedance  2(a).  The  resulting  modal 
equation  is  given  in  [3]  and  [4]: 


where 


M(a)  - MQ(a)  ♦ 4wtoZ(a)/k‘  =0  (1) 

Mo(°0  *=  -H^UjH)]  ♦P(S;H)  -a2Q(a;H)  (2) 


2 

r°°  exp(-UjH) 

j u ♦ u~ 

op 

P(a;H)= 

i? 

dX 

(5) 

2 

r°°  exp(-u^H) 

“2  2 

Q(c»;H)  = 

in 

dX 

(4) 

- n2Ul  +nlU2 

u • 

1 

2 2 4 

(\  - cj)  ; 

-n/2  < arg 

U1  < 

n/2 

(5a) 

U2  = 

2 2 4 

(x  ; 

-n/2  <.  arg 

U2  < 

n/2 

(5b) 

C1  * 

, 2 2,4 

(n  j - a ) ; 

0 < arg 

< 1\ 

(6a) 

^2  * 

, 2 2,4 

(n,  - a ) ; 

0 < arg 

< n 

(6b) 

o o 


Here,  Hv  J (x)  is  the  Hankel  function  of  the  first  kind  and  of  order 
o 

zero.  We  note  that  the  arguments  2 are  chosen  so  that  the  fields 
found  in  the  derivation  of  (1)  are  bounded  at  infinity  in  the  transverse 
direction.  The  arguments  of  u^  ^ are  chosen  so  that  all  of  the  integrals 
converge  everywhere  in  the  complex  X-plane.  While  equation  (1)  is  appli- 
cable to  a class  of  thin-wire  structures ^ ^ ^ we  are  specifically 
interested  in  the  application  to  a leaky  coaxial  line  composed  of  a center 

conductor  of  radius  c and  refractive  index  n ; a dielectric  insulator 

w 

around  the  conductor  of  index  n^;  a thin  uniform  metal  braided  sheath  of 

radius  b and  transfer  impedance  Z^,;  and  a dielectric  coating  around  the 

sheath  of  index  n and  radius  a as  shown  in  Fig.  2.  Provided  the  amount 
a 

of  leakage  is  small  , the  surface  impedance  can  be  found  by  calculating  the 
quantity 

Hlp.a) 

Z(c)  = — (7) 

2iTpH  (p,a) 

9 p=  a 

when  the  cable  is  driven  by  an  axially  symmetric  source,  where  is 

the  z-directed  electric  field  and  H is  the  <t>-directed  magnetic  field. 
According  to  Wait  [5],  Casey  [6]  and  Wait  ahd  Hill  [7],  this  transfer 
impedance  is  given  by 

Z(a)  = Za(a)  ♦ ZT(a)[Zb(a)  ♦ Z. (a)]/[ZT(a)  ♦ Zb(a)  ♦ Z.(a)]  (8) 


where 


Za(«) 


Zb(a) 


Z.(a)  = 


■>  2 2 

(n“  - a )k 
a o , ,a, 

2nTe_ uin7" 
o a 

, 2 2m,2 

(nb  )ko  , ,b. 

2nic  wii?  1 c 
o b 


it  k 

w o 

2nce  u>n2 
o w 


J(Uc) 

0 wo 

J . (t  k c) 

1 w o 
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The  equivalent  circuit  representation  of  (8)  is  also  shown  in  Fig.  2. 

We  note  that  Z . in  the  case  represents  respectively  the  series  impedance 

due  to  the  dielectric  layer  between  a and  b,  and  between  b and  c; 

Z^  is  the  internal  impedance  of  the  inner  conductor,  and  Z^  is  the 

shunt  transfer  impedance  of  the  braided  sheath.  For  a highly  conducting 

inner  conductor  Z.(a)  reduces  to 
l 


7 

* 

1 


ik 

o 


2iTce  ojn 

o w 


J (n  k c) 

0 w o 

J (n  k c) 

1 w o J 
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which  is  then  independent  of  a . 


The  transfer  impedance  of  the  braided  sheath  is  given  in  [5]  to  be 

2. 


ZT(a)  = 


-iko4 

e up 
o o 


(1  - 


a 


2 2 

n + n, 
a b 


) 


(12) 


where  is  the  transfer  inductance.  It  is  particularly  noteworthy  that 

the  total  surface  impedance  as  defined  in  (8),  is  a function  of  the 
propagation  constant  a . Furthermore,  it  is  easy  to  show  that  denominator 
of  (8)  may  be  rewritten  as 


where 


(ZT(a)  ♦ Zb(o)  ♦ Z.) 


Z'tcT 


“z5 


Z' 


ik  r . 

_L 

2ar.  u>  [_  2 
P nb 


4 


, 2 2. 
4(na  * V 


-] 


(13) 


(14) 


nb^na  * ntP  [v«0*-n (b/c)  + Lt  ♦ i2irZi/u>] 
t(na  * nb)v'olln(b/c)  + "bS1 


(15) 
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This  shows  that  the  Z(a)  has  a pair  of  poles  at  a = ±a^  . As  will 
be  shown  in  a later  section  of  this  report,  these  poles  have  a strong 
influence  on  the  location  of  one  of  the  roots  of  the  modal  equation. 

3.  Approximate  Expression  for  the  Integrals  P(ot;H)  and  Q(ot;H) 

The  integrals  P(a;H)  and  Q(a;H)  as  given  in  (3)  and  (4)  may  be 
computed  numerically  in  order  to  find  the  roots  of  M(ct)  . This  can  be  a 
time  consuming  process  considering  that  M(a)  must  be  evaluated  several 
times  for  each  root  found.  It  is  desirable  to  find  approximations  to 
P(a;H)  and  Q(a;H)  that  are  valid  in  the  regions  of  interest  and  are  more 
efficient  to  compute  than  direct  numerical  integration.  These  approxima- 
tions may  also  be  used  to  find  limiting  forms  of  M(a)  in  special  cases. 

In  [8],  approximations  to  P(a;H)  and  a slightly  different  form  of 
Q(a;H)  are  found  that  are  valid  in  the  region  |Cj|  <<  \^^\ . This  corre- 
sponds to  a in  the  neighborhood  of  n^ . Any  root  in  this  region  would 
be  highly  attenuated,  so  even  though  roots  may  exist  in  this  region,  they 
are  not  extremely  important.  Instead,  we  attempt,  in  Appendix  A,  to  find 
approximations  that  are  valid  in  the  region  |Cj|  » l^l*  because  any 
mode  found  in  this  region  has  a relatively  low  attenuation.  The  approxi- 
mations need  only  be  valid  for  | n^ ) H < 1 because  the  cable  usually  is 
placed  within  a skin  depth  of  the  surface  in  order  to  insure  sufficient 
penetration  of  the  wave  into  the  air  region.  In  such  a situation  we 
have  shown  in  Appendix  A that  P(a;H)  can  be  approximated  by 

P(a;IO  = j [4(1  -H.H)  ♦ (r  c,H/2N)24n(C,/t ;.)1  cxp(irH) 

iirUjHr  1 1 1 1 

* 2H<1)(C1H)  - 


(16) 
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where  is  the  Hankel  function  of  the  first  kind  and  order  m and 

ID 

only  the  principle  value  is  chosen  for  the  logarithmic  function. 
Approximation  for  Q(a;H)  is  somewhat  more  complicated  because  the 
integrand  as  given  in  (4)  has  a pair  of  poles  located  at  X « ±Xp  where 


Xp  - (aB  - o ) 


* aB 


,2  2,-»s 

n!n2(nl  + n2^ 


(17) 


and  1m  X^  > 0 in  the  complex  a-plane.  As  shown  in  Appendix  A,  the 
discontinuity  of  the  residue  contribution  at  X = X^  when  Xp  crosses  the 
real  axis  produces  a pair  of  square-root  branch  cuts  with  branch  points 
located  at  a ■ ia_  in  the  complex  a-plane.  An  approximation  which  specifically 
takes  this  singularity  into  account  is  obtained  from  (A. 41)  in  Appendix  A 
as 


2 2 

r 4i  nin2  2 A~l  2 2**  ' 

Q(a;H)  -I  - — jp — j exp (in^H/n)  / (a  -aRl  ♦ Q (a;ll) 

L n(n,  - n.)  1 


(18) 


Q'(a;H)  = -n^QlfajH)  ♦ n?<Ma;H)r  njQj(a;H)  (19) 

iTT(ni-n2)L 

where  n = (n^  ♦ n^)  . Expressions  for  Qj  for  j *»  1,2,3  are  derived  in 
(A. 27) , (A. 36) , (A. 39)respectively  in  Appendix  A and  arerepeated  here, 

Qj  (a;H)  = ittH^U.H)  - ( ~2$f  (exp(-ulpH)H<J 1 } (*iulj/C1  .C^O 
- exp(u  11)  Ho^  (iu^/Cj , CjH)] 


♦ PU  /X  ) cosh(Uj^H) (fcn(iUjp-iXp)  - tnCj] ; (20) 

Q2(a;ll)  *=  oxp(UjpH)Ein(ll[-iC|  -u^])  ♦ cxp(ulpH)l  in l»l [ - ic j ♦Ujp]) 

-[exp(pulpM)  *exp(ulpH)][ln(H)  ♦>]  ♦ texp(-ulpll)  - expiu^H)  «n(- iC^u  ) 

(21) 
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Qj(a;H)  = 2 exp(-ulpH){£n  2 -£n(-iC2)  - (u2p/Ap)  [An (-iu2p-iXp)  -An(-C2)]) 

2 2 i 2 2 i m (22) 

where  u,  = (A  -C,)  and  u = (A  -C  ) ; 1 is  the  Schwarz 

lp  p 1 2p  p Z x 

or  Lipschitz-Hankel  integral  defined  in  [12]  as 


(a,z) 


(1) 

exp(iat)HUJ  (t)dt 

'n 


and  can  be  computed  from  a series  expansion  given  in  (A. 23).  Ein  is  the  modified 
exponential  integral  given  in  [10]  as 


Ein  (z) 


k = l k(k!) 


Finally,  the  logarithmic  function  An(z)  is  specified  by  its  principle 

value.  It  should  be  noted  that  the  behavior  of  Qr(a;H)  given  in  (19) 

is  not  singular  at  a = , and  its  value  at  may  be  found  using 

(A. 27b) , (A. 36b)  and  (A. 40b).  However,  the  first  term  in  (18)  does  contain 

2 2 

explicitly  the  term  (a  -ctg)  so  that  the  expression  of  Q(ot;H)  will 
blow  up  when  a = c<g. 


4 . Classification  of  Modes 

Because  the  leakage  of  a braided  cable  is  usually  very  small  (i.e. 
L.j.«uo/k  ) it  is  logical  to  view  the  propagating  modes  of  a buried  leaky 
cable  as  resulting  from  coupling  of  modes  between  the  coaxial  cable  and  the 
external  waveguiding  system  of  a buried  and  insulated  wire  located  near 
the  earth  surface.  Thus,  by  inserting  (8)  and  (13)  into  (1),  t^e  modal 
equation  can  be  rearranged  to  reflect  such  a viewpoint: 

(ot2  - a ^)  [Mq  (a)  ♦ (4a>eo/k2)Za  (a)  ] + A (a)  = 0;  (23) 


r 


n 


whore 


A (a)  - (^)(zb(a)  ♦ zi(aJJzT(«)/z' 


(24) 


and  |A(a)|  «1  because  Z is  assumed  to  be  small.  For  simplicity. 


we  let  Z^(a)  “ 0 so  that  the  center  conductor  of  the  leaky  coaxial  is 


perfectly  conducting.  As  is  evident  in  (23)  and  the  definition  of  Zy 


in  (12),  if  we  now  let  the  transfer  inductance  of  the  braid  E approach 
zero  then  we  have  o * c»7  » n^,  which  obviously  corresponds  to  the  TEM- 
mode  of  a coaxial  line  with  inner  and  outer  radii  c and  b.  Provided 
the  leakage  remains  small,  we  can  find  the  solution  of  the  modal 
equation  by  perturbation;  the  zeroth-order  solution  is  ot  ■ ct^-;  the 
first -order  solution  is  found  by  inserting  a , for  a in  M (a),  Z (a) 

“ Oft 

and  A(oi);  the  second-order  solution  is  found  by  inserting  the  first- 


order  solution  in  M (a),  Z (a)  and  A(ot);  etc.  Similar  to  the  TEM-mode 

O ft 


in  a coaxial  line,  we  expect  the  current  in  this  case  to  be  almost  equal 
and  opposite  on  the  inner  conductor  c and  the  inner  surface  of 
conductor  b.  Such  a mode  is  designated  as  a bifilar  mode.  On  the  other 
hand,  other  acceptable  solutions  of  (23)  when  A(a)  ► 0 may  be  found  from 


M0(ot)  ♦ (4u*co/k^)Za(m) 


0 


(25) 


which  obviously  represents  the  modes  supported  by  an  insulated  conducting 
wire,  or  a Goubau  line  buried  near  the  earth  surface,  and  with  a surface 
impedance  of  Z . As  shown  in  [2],  such  a line  supports  at  least  one 
mode  with  a known  propagation  constant  in  the  low  frequency  limit.  The 
current  is  now  highly  concentrated  on  t he  outer  surface  of  conductor  b, 
and  consequently,  the  mode  is  designated  as  a monofilar  mode.  Obviously, 
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tho  value  of  a may  be  found  again  from  a perturbative  scheme  based 
upon  a zero-order  solution  obtained  from  (25). 

In  addition  to  the  monofilar  mode,  it  is  known  that  a second  mode 
near  o = also  exists  in  the  case  of  an  elevated  Goubau-line  in  air 
and  located  near  the  earth  surface  [9],  [12).  Mathematically,  the 
occurrence  of  this  mode  is  heavily  influenced  by  the  inverse  square- root 
singularity  of  the  Q- integral  in  the  modal  equation;  such  a singularity 
is  displayed  explicitly  in  the  approximate  form  of  Q given  in  (18). 

The  dominance  of  the  Q-integral  also  means  physically  the  field  distribu- 
tion of  such  a mode  will  be  much  more  spread  out  along  the  air- earth 
interface  than  that  of  the  monofilar  mode,  and  consequently,  is  designated 
as  a surface-attached  mode.  In  Appendix  B,  it  is  shown  that  this  mode 
indeed  can  exist  for  a sufficiently  small  buried  depth  H,  but  may  cross 

2 2 i 

the  branch  cut  (Re (a  - a^)  =0)  and  become  an  improper  mode  for  a 

larger  H. 

Before  we  present  the  numerical  evaluation  of  modes  for  a general 
case,  another  special  case  of  interest  is  when  the  wire  is  located  at 
the  air-earth  interface.  Consider  tho  limit  as  H approaches  zero, 
while  keeping  the  ratio  A:H  sufficiently  small  so  that  the  thin  wire 
approximation  continues  to  hold.  If  we  keep  the  relative  cable  geometry 
the  same  then  the  transfer  impedances  Z (a),  Z (a) , Z (a) , and 
quantity  Z'  will  remain  constant.  We  will  again  assume  that 
Z^(a)  is  zero  for  simplicity.  The  two  llankel  functions  in  the 
expression  for  Mq(oi)  in  (2)  diverge  individually;  however,  their 
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difference  approaches  a constant,  By  taking  the  limit  as  H approaches 
zero  in  (16)  and  (18)  we  find  that  both  P(a;H)  and  Q(n;H)  diverge 
logarithmically,  hence  in  this  limit  MQ(a)  becomes 


V“’  ' - °2J 


Tifnj  ♦ n2) 

By  inserting  (26)  into  (23)  we  find  that  we  can  neglect  Z (a),  Hence 


(26) 


we  obtain 


2 2 2 2 
_ „ /n  riTi(n,  ♦ n-)i 

(c2  - c,  )(.-4^  - A = — s— i.  | a (a) 

x ' 1 4 In  H J 


(27) 


As  H approaches  zero  the  term  on  the  right  side  of  (27)  vanishes  so  in 
this  limit,  the  value  of  a for  the  monofilar  mode  reduces  to 

a = [ (n^  ♦ n2)/2]^ 


which  agrees  with  the  well-known  result  obtained  by  Coleman  [11]  for  a 
thin-wire  located  in  the  air-earth  interface. 


5.  Numerical  Results 

We  have  developed  a computer  program  to  compute  the  roots  of  the 
modal  equation  (1).  This  program  computes  P(cx;H)  and  Q(a;H)  either  by 
direct  numerical  integration  of  (3)  and  (4)  or  by  the  approximations 
given  in  section  3.  Unless  specified  otherwise,  the  refractive  index 
of  the  earth  is  taken  to  be  = S.3  ♦ i0.95,  and  for  the  air  n2  = 1.0. 
The  frequency  is  f = 100  MHz,  and  the  dimensions  of  the  cable  in 
reference  to  figure  2 are:  a = 1.15  cm;  b = 1.0  cm;  c = 0.4  cm.  The 
braid  inductance  is  Lj,  * 40  nil  , and  the  inner  conductor  is  assumed  to 
be  perfectly  conducting  so  that  Z^  = 0.  The  refractive  index  of  the 


r:  ..  . 
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• i 

coating  was  n a 1.449  which  corresponds  to  Teflon.  The  insulator 
a 

refractive  index  was  varied  from  n^  = 1.0  to  1.449,  and  the  depth  at  which 
the  cable  was  buried  was  varied  from  h=>0,l  m to  1.0  m.  The  skin  depth 
of  the  earth  is  about  0.5  m so  that  the  root  locations  should  be 
relatively  independent  of  h when  h > 0.5  m.  In  all  cases,  we  found 
two  and  sometimes  three  distinct  roots  for  any  given  set  of  parameters. 

The  first  mode  is  a bifilar  mode  in  which  the  currents  on  the  inner 
conductor  and  the  braid  are  approximately  equal  but  have  opposite  signs. 

This  mode  has  the  least  attenuation  of  any  of  the  modes  since  the  fields 
are  concentrated  on  the  inside  of  the  cable.  The  location  of  this  root 
tended  to  follow  cx^  when  the  value  of  n^  was  changed.  A plot  of  the 
location  of  this  mode  as  a function  of  h is  given  in  figure  3. 

The  second  mode  is  the  monofilar  mode  in  which  almost  all  the  current 
is  on  the  braid,  so  this  mode  is  similar  to  the  case  of  a buried  coated 
wire.  A plot  of  the  location  of  this  mode  is  given  in  figure  4.  Note 
that  this  mode  has  an  attenuation  of  about  12  to  15  dB/m.  Although  it 
is  heavily  attenuated  along  the  line,  such  a mode  has  the  major  part  of 
field  distribution  located  outside  of  the  cable  and,  hence,  is  capable 
of  interacting  with  surrounding  objects  in  earth. 

The  third  mode  is  the  surface  attached  mode  which  has  its  fields 
concentrated  near  the  interface  between  the  air  and  earth.  As  shown  in 
figure  5 this  mode  only  exists  for  certain  values  of  h.  For  h > 0.34  or 
h < 0.19  m,  the  root  crosses  the  branch  cut  onto  the  improper  Riemann  surface 
in  the  complex  a plane.  The  mode  becomes  improper  in  these  cases  and  is 
absorbed  in  the  surface  wave  radiation.  Even  when  this  mode  does  exist 
it  is  located  very  close  to  the  branch  point  a given  by  (17)  . For 

h 


.00572 


1.449,  ru  = 1.449,  L_  = 40  nH,  f = 100  MHz. 
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this  reason  the  fields  of  this  mode  will  b.e  very  similar  to  the  ground 
wave  fields  which  are  spread  out  along  the  entire  surface. 

The  accuracy  of  the  approximations  to  P(ot;H)  and  Q(a;H)  was  detexmined 
by  comparing  the  values  of  the  roots  found  using  the  approximations  to  the 
values  obtained  using  the  numerical  integration.  For  the  bifilar  mode 
at  h = 0.1  n the  values  of  the  roots  agrees  to  within  10  J and  at  h=0.1m, 
they  agreed  to  within  8 * 10  **.  For  the  surface  attached  mode  the  roots 

_9 

agreed  within  10  for  all  values  of  h.  For  the  monofilar  mode  the 

approximations  gave  completely  inaccurate  results  except  for  very  small  h 

say  h <.01  m,  because  the  assumption  that  |c^|  >> | ^2 1 no  l°nger  valid 

in  the  region  where  a of  this  mode  is  located.  Thus,  only  direct  numerical 
integration  is  used  in  this  case. 

6.  Concluding  Remarks 

In  this  report  we  have  investigated  the  roots  of  the  modal  equation 
for  a leaky  coaxial  cable  buried  in  a lossy  earth.  We  have  found  that 
three  distinct  modes  do  exist  for  most  cable  depths  of  practical  interest. 
We  have  given  perturbation  formulas  to  locate  the  roots  of  two  of  these 
modes . 

Of  the  three  modes  the  bifilar  mode  is  probably  the  easiest  mode  to 
excite  since  the  fields  are  similar  to  a TEM  mode  in  a lossless  coaxial 
cable.  This  mode  has  the  least  attenuation  but  also  has  only  a small 
amount  of  field  penetration  into  the  air  region  because  the  total  current 
on  the  cable  is  nearly  zero.  The  surface-attached  mode  has  more  field 
strength  in  the  air  for  a given  amount  of  current  on  the  inner  conductor 
than  the  bifilar  mode.  For  this  reason,  such  a mode  is  ideally  suited  for 
the  design  of  a groundwave  detect i on  system.  But  since  the  location  of 
this  mode  is  close  to  the  branch  point  at  a = one  may  not  be  able  to 


.! 
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excite  such  a mode  independently  without  substantial  surface  radiation. 

The  monofilar  mode  has  fields  concentrated  in  the  earth  region  so  that 
this  mode  is  highly  attenuated. 

From  the  above  discussion  it  is  quite  obvious  that  there  are  several 
conflicting  considerations  in  the  design  of  a wave  guiding  system  for 
detecting  anomalies  (or  intruders).  On  the  one  hand,  the  bifilar  mode  is 
easy  to  excite  and  once  excited,  it  is  capable  of  propagating  along  the 
line  with  least  amount  of  attenuation.  However,  its  field  distribution 
is  highly  concentrated  to  the  interior  coaxial  region  so  that  it  is  not 
particularly  sensitive  to  anomalies  located  above  the  earth  surface.  On 
the  other  hand,  the  monofilar  mode  suffers  higher  attenuation  and  the 
surface-attached  mode  can  not  be  excited  easily.  The  question  then  is 
how  to  bring  about  an  optimized  system  which  provides  the  best  compromise 
to  these  conflicting  considerations.  To  this  end,  it  may  be  useful  to 
first  define  a performance  index  of  a specific  mode  as 

P.I.  = £ (28) 

PTIm(a) 

where  .°° 

6 = j |E  | ^dy  is  the  integration  of  the  magnitude  square 

^oo  n 

of  the  vertical  electric  field  component  along  the  earth 
surface  at  a given  cross-section; 

P.p  is  the  total  amount  of  time-average  power  flow  across  the 
same  cross-section; 

Im(a)  is  the  imaginary  part  of  a or  the  attenuation  constant. 


" 


Although  such  a definition  is  not  at  all  unique,  it  is  reasonable  to  assume 
that  the  set  of  parameters  that  yields  the  highest  performance  index 
according  to  (28)  provides  the  "optimum"  design  for  a detection  system. 
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In  Appendix  C expression  for  the  field  components  for  each  discrete 
mode  are  given.  However,  the  actual  computation  and  intercomparison  of 
the  performance  indices  of  various  modes  will  be  presented  separately  in 
a supplemental  report  at  a later  date. 
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Appendix  A 

In  this  Appendix,  approximate  expressions  are  derived  for  the 
integrals  P(a;H)  and  Q(a;H)  under  the  assumption  that  |c  | » | C2 1 and 
| n^ | H < | n j | H <1.  It  is  easy  to  see  from  (3)  that  the  main  contribution 
to  the  P-integral  comes  from  the  range  of  X where  u^  can  be  approxi- 
mated simply  by  X . We  therefore  divide  the  integral  in  (3)  into  two 
parts , 

P(o;H)  = Pj  Cci;H)  ♦ P2(a;H)  (a.1) 


where  P (a:H)  represents  the  dominant  contribution,  i,e,, 

4 exp(-u^H) 


r expv-c 
>,  (a;H)  = 

1 111  J u + 


dX 


o 1 


r r -UjH 

Xexp(-UjH)dX  - u^e  dX 


l^2l  10 


which  is  known  analytically  to  be 


Pj (a;H)  = 


iTTC, 


H „2 


exp(i^jH) 


*2II01)(C1H)  - H1(1)^1H) 


(A. 2) 


The  remainder  term  P2(ci;H)  is  now  given  by 

00 

p2 H f [u^2  - S^TxJ  exp(-UlH)«lX 


(A. 3) 
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Note  that  the  term  in  the  brackets  decays  as  X'  for  large  X so  the 
important  part  of  the  integral  where  X is  small.  For  H small  enough 
the  exj  "■>■£  »;ial  term  is  constant  for  X less  than  some  Xq  so  the  exponen- 
tial term  may  be  replaced  by  its  value  at  X = 0 for  the  integration  of  X 
between  0 and  Xq.  Since  the  integrand  decays  as  X ^ for  large  X , 
the  integration  from  Xq  to  » is  insignificant  even  without  the  exponen- 
tial decay.  This  allows  us  to  let  approach  infinity  without  affecting 
the  results,  The  exponential  term  is  replaced  by  its  value  at  X = 0 
instead  of  unity  so  that  P2(a;H)  remains  small  when  compared  to  P (« ; H ) 
as  n^  approaches  infinity.  Hence  P2(a;H)  is  approximately 


XQ 

>iH)  = H exP(iciH)  “1  ( [57VSJ-  uj— x]  dX 


expCi^HHn^/Cj) 


(A. 4) 


2 2 2 

where  N = nj  " n2  » the  principal  branch  of  the  logarithm  is  chosen. 

Substituting  (A. 2)  and  (A. 4)  into  (A.l)  yields 


P(a;H) 


r ;2alH  ( C2^ 

T 1 -UjH  + 5 — *n  (— p— ] exp(ic.H) 

H L 1 2NZ  V W 1 


(A. 5) 


This  approximation  for  P(a;H)  is  valid  for  |c2l  <<1^1  or  for  H|^j|  <<  1. 
In  the  latter  case  (A. 5)  may  be  further  approximated  by  using  the  small 


argument  expansion  of  the  Hankel  functions  to  obtain 

2i  i ^2*  (**2  \ 

P(a;H)  = ^ [*n(C,H)  + y - 2]  + 1 - ± In  ' f- ) 

TT  i " TTN^  '^1  ' 


(A.  6) 
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where  y = 0.577216  is  Euler's  constant.  The  result  Riven  in  (A. 6) 

agrees  with  the  approximate  formula  derived  earlier  by  Chang  and  Wait  in  [2]. 

On  the  other  hand,  by  multiplying  and  dividing  the  integrand  by  the 
2 2 

factor  OjU^  - n2Uj»  tho  integral  Q(a;H)  as  given  in  (4)  may  be  rewritten 
as 


Q(a;H) 


itr(n 


2 f ^ 

Sni)  J ' 


2 2 
n.u2  - u2ul 


a2‘ap 


exp(-UjH)dA 


(A. 7) 


where 


2 2 ^ 

Xp  « CaB  - a ) 

- (Cj  “ n*/n2)* 

n ■ (n~  ♦ n2)* 

aB  * nlV" 


0 i arg  X < v 
P 

0 <„  arg  n < a 


(A.  8) 
(A.  9) 
(A. 10) 


The  integrand  in  (A. 7)  appears  to  have  a pair  of  poles  at  X = ±A^,  the 

location  of  which  are  a function  of  a.  According  to  [8]  the  discontinuity 

of  tho  residue  calculation  at  X ■=  X when  X crosses  the  real  axis  causes 

P P 

a branch  cut  in  the  ci-plane  of  Q(a;ll).  In  addition,  Q(oi;H)  becomes 

unbounded  as  a approaches  the  branch  point  because  X ^ tends  to  zero 

and  the  integral  in  (22)  does  not  converge  when  X^  equals  zero. 

In  order  to  insure  that  tho  poles  do  exist,  the  numerator  in  (A. 7) 

needs  to  be  evaluated  at  X * X , and  verified  that  it  is  non-zero  there. 

P 

In  the  physical  case  that  region  1 is  the  earth  and  region  2 is  air, 
tho  following  relationships  hold: 


arg  iij  > arg  n > arg  n^  > 0 


(A. 11) 


Under  those  constraints  it  can  be  shown  that 


J 
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2 * 

uil  B ui„  “ -in./n 
^X-X  XP  1 


(A. 12) 


If  in  addition  we  can  assume  that 

A 

arg  n > 2 arg  n~ 


(A. 13) 


then  it  can  also  be  shown  that 

2 

U2|  = U2p  “ in2'n  CA. 14) 

•x=x 

P 

Using  (A. 12)  and  (A. 14)  the  numerator  of  the  integrand  of  (22)  is 

evaluated  at  A * A to  be 
P 


2 2 
nlU2  “ n2Ul 


7 2 A 
= 2 in^n2/n 


(A. 15) 


X^A 


which  is  non- zero,  hence  the  integrand  does  have  poles  at  X = ±X^.  Note 
however  that  the  relationship  (A. 13)  may  not  always  hold  if  both  regions 
1 and  2 are  lossy,  in  which  case  the  sign  in  (A. 14)  should  be  revorsed 
and  the  poles  would  disappear. 

An  approximation  for  Q(a;H)  will  now  be  found  in  which  the  singu- 
larity at  a * a is  accounted  for  exactly.  In  addition  the  approximations 
will  have  the  proper  limit  as  approaches  zero,  or  as  11  approaches 
zero.  Rewrite  (A. 7)  as 


Q(a;H) 


where 


-f—j  {-  n^Q(a;H)  + njQ2(a;H)  + njQ,(a;H)) 


in  (n^  - n2) 


a> 

Q^ot;!!)  - j -^y-expt-UjHJdX 

* **p 

Q2 (a;  11)  - 2 J • expf-u^OdX 


X-X, 

o P 


(A. 16) 


(A. 17) 


(A. 18) 
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u„  - A 


Q3(a;H)  = 2 2 exp(-u^H)dA 


A -A 

CO  p 


The  integral  QjCajH)  is  evaluated  as  follows. 


a2  ir*  expC-u  H) 
Qj(cilH)  = ^ ) I 2 2 

3H  I*-  (ul  ‘ ulp)ul 


dA 


2u  2 ^11  ^a’^’ulp^ 

/ulp  9fT  11  ip 


Qll(a;h'-Ulp)} 


where 


C exp(-UjH) 

Qji  (a;H»uip)  - j (uru  )uj 


dA 


IP 


[expC-CUj-u  )H)  - 1] 


«iexpi-iu  -u 
7 LTi£ 

(uruiP)ui 


dA 


w 

+ f _dA__ 

.1  <vVui . 

j t°  exp(-u  H) 

exp(-uipH)  j -jo expCiiiptjj  — - 

' I 


dA  dt 


dA 


A ‘-A 

-OO  p 


(a'-A^i 


= exp( 


-vfr 


,(D 


iff 


exp(uipt)Ho  '(CjtJdt  ♦ j- 


U1P  \ 


dA 


(A2- A2)u  I 

-CO  V J J . 


(A. 19) 


(A.  20) 


(A. 21) 


In  order  to  evaluate  the  finite  integral  in  (A. 21)  we  can  use  the 
incomplete  Lipschitz-Hankel  integral  or  Schwarz  function  which  is  defined  as 


I 
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ittu 


QjfojH)  = Qj(ct;H)  ♦ -y^  exp(-ulpH) 


(A. 26) 


where 


ITTU 


1U, 


Qj(a;H)  - - exp(-ulpH)He(1)(-  , ^H) 


ITTU. 


♦ ~2~^  exp(ulpH)He(1)(-^£,  ^H)  ♦ iwH^KjH) 


+ [exp(-ulpH)  ♦ exp(ulpH)]  [An(iulp  - iAp)  - «.n  Cj ] 


IP 


lp  P 


(A. 26a) 


The  second  term  on  the  right  side  of  (A. 26)  contains  all  of  the 
singular  parts  of  Qjta.H)  near  the  branch  point  afi.  At  a = aB , Q^o.H) 


is  finite  and  reduces  to 


Q|  (°b»H)  - iirH^^iu  .H)  ♦ uipHllHi(1)(iuipH) 


lp  ' lp 


IP 


(A. 27) 


The  integral  Q2(a;H)  may  be  evaluated  as  follows: 
00 

Q2(a;H)  = 2 f -J-— - exp(-u  H)dA 
J u,  -u. 


0 1 lp 

1 


u,_  ^22^a’H,Ulp^  " ^22^a;H,"Ulp^ 


(A. 28) 


IP 


where 


oo 

Q22Co;H.ulp)  ■ J oxpt-UjWdx 

0 p 

oo  00 

= exp(-ulpH)j  exp(ulpt)  j X expC-UjtJdA  dt  (A. 29) 


H 0 

This  last  step  is  valid  only  if  Real(u1~u1  ) > 0 for  all  X . This  is  true 


if  Real(-iC1  -ulp)  > We  wil1  derive  an  expression  for  Q22(tt;H,ulp) 


IP' 


1 
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assuming  that  a is  in  the  region  in  which  the  above  inequalities  hold. 
We  will  then  analytically  continue  the  resulting  expression  for  all  values 
of  a . Note  that  this  problem  does  not  exist  in  the  derivation  of 
(^(oijH.-Ujp)  because  Real(u^  + Ujp)  > 0 for  all  X and  aH  a • 

Doing  the  X integration  in  (A. 29)  yields 

00  r I 

Q22(«;H,u  ) - exp(-ulfH)j  L-r^  + “H  ®xp[(iC1  + ulp)t]dt 

LI  t 


- jjexpCi^H)  + u exp(-u  H)E1(H[-iC1  -u  ]) 


(A.  30) 


Ej(Z)  and  E^CZ)  are  the  exponential  integrals  of  order  1 and  2 respectively 


[9].  Alternately  E^(Z)  may  be  expressed  as 


Ej (Z)  = Ein(Z)  - A n(Z)  - y 


(A. 31) 


where  y is  Euler's  constant  (=0.577216),  and  E^(Z)  is  an  entire 


function  which  has  the  expansion  [9] 


Ein(Z) 


k-,1  k k! 


(A. 32) 


Inserting  (A. 31)  into  (A. 30)  we  obtain  an  expression  for 


Q22(a;H,ulp) 


Q22(a;H,ulp)  * pf  expCi^jH)  + ulp  exp(-uJpH)Ein(H[-ic1  -ulp]) 
- ulp  exp(-u  H)[In  H+fcnf-i^  - u ) + y] 


= i expCiCjH)  ♦ ulp  exp(-ulpH)Din (H[-iC^- ulp] ) 

- ulp  exPC-ulpH)[tn  II-  ln(-iCj  *ulp)  ♦ y] 

’ 2 ulpexp('ulpH)*n(_iV 


(A. 33) 


(A. 34) 
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The  only  singularities  in  (A. 34)  are  the  branch  cuts  due  to  the 
square  root  in  Cj  and  the  logarithm  of  Xp,  hence  this  form  is  the  proper 
analytic  continuation  for  all  a. 

Inserting  (A.  34)  and  (A.  33)  with  -u^  substituted  for  u^  into  (43) 
we  obtain  the  following  expression  for  Q2(a;H) 


Q2(a;H)  = Q2  (a ;H)  - 2 exp(-ulpH)Jln(-iAp) 


(A. 35) 


where 


Q2(a;H) 


eXp(_UlpH)EinWiCl  'ulpJ)  + expCu^HjEinCHf-iCj  +ujp]) 


- [exp(-ulpH)  ♦ exp(ulpH)]  [fcn(H)  + y] 

♦ [exp(-ulpH)  - expCu^HJjant-iCj  + ulp)  (A. 36a) 


Q2(a;H)  has  a logarithmic  singularity  at  a = a R due  to  the  last  term 
in  (A. 35).  Q2(a;H)  has  no  singularity  at  a = afi  and  may  be  evaluated  at 

this  point  to  obtain 

Q2(aB;H)  - exp(ulpH)Ein(2ulpH)  - (exp(-ulpH)  ♦ exp(ulpH) ] [*n(H)  +Y] 

♦ [exp(-ulpH)  - exp(ulpH)]fn(2uJp)  (A. 36b) 

The  integral  Q3(a;H)  may  be  rewritten  here  as 

j [°  exp(-u  H) 

Q (a;H)  - -2C  dA  (A. 37) 

2 Jo  (u2*A)(A2-A2) 

From  this  expression  it  is  apparent  that  Q3(a;H)  vanishes  as  C2  approaches 
zero.  The  integral  also  converges  even  when  H is  set  to  zero  whereas  both 
QjCajH)  and  Q2(a;H)  diverge  logarithmically  as  H approaches  zero. 
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We  therefore  expect  the  major  part  of  the  integration  to  come  from  the 
region  near  the  poles  ±Xp  in  the  complex  X-plane.  Since  the  integrand 
decays  as  X for  large  X , we  will  approximate  Q3(o;H)  by  replacing  the 
exponential  term  by  its  value  at  X = X . Hence  Q,(oi;H)  becomes 


o 

Q3(a;H)  = 2 exp(-ulpH)  j 


U2  X 

2 2 
X -X 


dA 


= 2 


V 

fX0 

oo 

l &- 

f XdA  . .2 

f dX 

ll  “2 

o x2-x2  U2P 

l Mx2-xf)J 

= 2 exp(-ulpH)js.n(2)  - «,n(-ic2)  ♦ fcnt-iAp) 


t 

2 

u_ 


dX 


(A. 38) 


The  infinite  integral  in  (A. 38)  is  given  by  (A. 24)  with  C substituted 

r 2 2 * 

for  tj.  We  may  also  replace  (c2  - Xp)  with  -iu2p  so  that  Q3(c*;H) 
becomes 


itiu 


Q,(a;H)  = Q' (a;H)  + - 2 Jtn(-iX  ) > exp(-u.  H) 


IP 


(A. 39) 


where 


Qj(a;H)  = 2 exp(-ulpH)  j*n(2)  - Hn(-i<;2)  - ^[*n(-iu2p-iXp)-Jln(-i;2)  ]| 


(A. 40a) 


Q3(a;H)  does  not  have  a singularity  at  a = cig  and  may  be  evaluated  at 
this  point  to  obtain 


Q;(Vh)  = 2 cxp(-ulpH)  [ JLn (2)  -£n(u2p)  - 1] 


(A. 40b) 
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The  logarithmic  singularity  at  a = ag  in  the  last  term  of  (A. 39) 
will  exactly  cancel  the  singularity  of  (^(ajH)  in  (A. 35).  By  inserting 

(A. 26),  (A. 35)  and  (A. 39)  into  (A. 16)  we  may  write  Q(a;H)  as 

2 2 2 

4i  n^n  in*H 

Q(«;H) A 4--4-  exp(-^—  ) ♦ Q'(a;H)  (A.41) 

Xpn(nj-n2)  n 

where 

Q'(a;H)  = {-n^Q^ajH)  + n^(a;H)  ♦ n^(a;H)}  (A.42) 

inC^-^) 

Q|,  Q2,  and  are  defined  in  (A. 27),  (A. 36)  and  (A. 39)  respectively. 

It  should  be  noted  that  q' (aJH)  is  not  singular  at  a = ag  and  its  value 
at  this  point  may  be  found  using  (A. 27),  (A. 36b),  and  (A. 40b). 
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Appendix  B 


In  order  to  demonstrate  that  a surface-attached  mode  may  exist  in 
the  buried  cable  problem,  we  first  insert  (18)  into  (1)  and  (2)  to  yield 
an  alternative  form  of  the  modal  equation  as 


4ia2n2n2  exp(in2H/n) 


p n(n* -n*)  [nT  (a)  + 4weoZ(a)/k2] 


(B.l) 


where 


M;ca) 


(Cl/ni)  I”®15  (^lA)  "Ho1)(5iH)l  + pk;H)  -cVftiW)  (B.2) 


and  Q' (a;H)  is  defined  in  (19).  The  right  side  of  (B.l)  is  a smooth 

function  of  a near  a , so  we  may  insert  a for  a and  calculate  X . 

B B P 

We  may  then  calculate  a using  the  expression 


2 2 . 2 2 2,.  2 2.  .2 

« = aB  ' Xp  = nln2/(nl  + n2}  ' Xp 


(B.3) 


This  value  of  a may  be  inserted  back  into  (B.l)  to  recalculate  a more 
accurate  value  for  a . 

The  imaginary  part  of  X must  be  positive,  so  that  if  the  right 

P 

side  of  (28)  has  a negative  imaginary  part  then  the  mode  does  not  exist. 

We  have  found  that  this  depends  on  the  value  of  H.  In  the  limit  as  H 

approaches  zero,  (B.l)  at  a = a reduces  to 

B 


_44 
2im^  n^ 


n(n^  -n2)(n^  + n2Hn(H) 


; n = (n2  ♦ n2)  . (B.4) 
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If  I n ^ | is  greater  than  |n2|  and  they  satisfy  the  condition  that 


arg  n > 2 argn  , then  the  right  side  of  (B.4)  has  a positive  imaginary 


part.  So  for  very  small  H a root  of  the  modal  equation  exists  that  is 


close  to  a = aR.  As  H is  increased,  this  root  may  or  may  not  exist 


depending  on  the  imaginary  part  of  the  right  side  of  (B.l).  In  all  cases. 


as  will  be  shown  in  a later  section,  when  the  root  does  exist  it  is  very 


close  to  a » cu. 
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Appendix  C 


In  this  Appendix,  expressions  for  the  fields  of  the  discrete  modes 
are  given.  To  be  consistent  with  the  thin  wire  approximation,  all 
higher  order  terms  of  kQa  are  dropped.  Generalizing  the  results  of  [ 3] 
we  can  express  the  fields  in  terms  of  the  z-components  of  the  electric 
and  magnetic  Hertz  potentials  to  obtain 


Ex(x.y»a)  = iko[a3U1>2/ax  + 2/3y] 


Ey(x»y»a)  = ikota3U1  2/3y  - n^Vj  2/3x] 
Ez(x,y,a)  = k*  U12 


Hx(x,y,a)  = iko[ct3V1>2/3x  - (n;f2/n0)8Ulf2 


Hy(x,y,a)  = ik0[a3V12/3y  + (n^  2/r>0)3Ui  2 


H,(x.y.a)  - <?li2  k2vJ2 


/3y] 

/3x] 


The  subscript  1 is  used  for  the  earth  region  (x  > 0)  and  2 is  for  the 
air  region  (x  < 0),  see  Figure  1.  The  potentials  are  found  in  [1]  or  [2] 
to  be 


-n  I 

U.  (x,y,a)  = — j2— 

4^Iko 


\r2 

2 

ff  "in[* 

:,k0(tx-h)2  * r2)  ] ‘ «'*>[ 

l"i  L 

"l 

2 2 ' 

r*y  ) 


* £ 


a 


♦u,  2 2 

2 niu2  ♦ n2Uj 


exp(-u1kQ(x+h)  - iXkQy)dX 
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U2(x,y,oQ  = 


V 


2uk  C_  m 
o 2 00 


LVU2 


1 


2 2 J e*PC-u1k0h  ♦u2kQx-iXkoy)dX 


nlU2  + n2Ul 


V^x.y.a) 


V2(x,y,a)  = 


al 

1 

ni 

X 

J 

00 

lVu2 

"lU2  + 

2 

n2Ul 

U1 

al 

00 

r 
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2 

ni 

1 

X 

2TtkoC22  J 

00 

L v«2 

2 2 
nlU2*n2Ul 

U1 ' 

In  these  expressions  I is  the  total  current  on  the  cable.  2 2 

are  defined  in  (5)  and  (6). 
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